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We consider random sequential adsorption processes where the 
initially empty sites of a graph are irreversibly occupied, in ran¬ 
dom order, either by monomers which block neighboring sites, or 
by dimers. We also consider a process where initially occupied sites 
annihilate their neighbors at random times. 

We verify that these processes are well defined on infinite graphs, 
and derive forward equations governing joint vacancy/occupation 
probabilities. Using these, we derive exact formulae for occupation 
probabilities and pair correlations in Bethe lattices. For the block¬ 
ing and annihilation processes we also prove positive correlations be¬ 
tween sites an even distance apart, and for blocking we derive rigorous 
lower bounds for the site occupation probability in lattices, includ¬ 
ing a lower bound of 1/3 for . We also give normal approximation 
results for the number of occupied sites in a large finite graph. 


1. Introduction. This article is concerned with certain Markov processes 
with state space {0,1}^, where V is the vertex set of a graph, in which the 
evolution at each site of V is monotonic in time and the stochastic rate 
of evolution at a site is either a constant or zero, depending on the local 
configuration. Unusually for interacting particle systems, the focus is not 
on equilibrium measures (owing to the irreversibility) but rather on exact 
computation or estimation of probabilities, in which more progress than 
usual can be made because of the simple form of the interactions. 

Two of the processes which we consider are versions of random sequential 
adsorption (RSA), in which the graph represents an initially “vacant” sur¬ 
face onto which particles are deposited successively at random, subject to 
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an excluded volume effect whereby a deposited particle prevents new parti¬ 
cles from being deposited nearby. We consider monomer RSA with nearest- 
neighbor exclusion/blocking (or blocking RSA for short) and dimer RSA. We 
also consider two forms of annihilation process in which sites are initially 
alive and from time to time kill their neighbors. 

Interest in RSA from the physical sciences community is considerable: 
see the 1993 survey of Evans [7], and more recently, the series of surveys 
in Volume 165 (2000) of the journal Colloids and Surfaces A. The principal 
motivation comes from the notion of particles being physically deposited 
onto a surface. Motivation is discussed further in [19], where qualitative 
results are given for a class of process which generalizes those considered 
here on the integer lattice. Continuum versions of RSA are also of interest 
(see, e.g., [18, 20]) but are not considered here. The annihilation process was 
introduced by O’Hely and Sudbury in [14], as a model for the thinning of 
seedlings, where a plant can grow only by eliminating its neighbors. 

In the present paper we consider general graphs (not only lattices). On 
the qualitative side, we derive an existence result for these processes on 
infinite graphs, forward equations governing the rate of change and a normal 
approximation theorem for the sum of the states of the vertices of a large 
finite graph. On the quantitative side, we derive explicit formulae for the 
occupation probability of sites in regular infinite trees (Bethe lattices) and 
also for correlations between sites in such trees. Moreover, we show that 
in the case of blocking RSA, these formulae provide good rigorous bounds 
for occupation probabilities in lattices snch as comparing favorably to 
previously obtained rigorons bounds. The proof of these results for blocking 
RSA uses a result on positive correlations between sites an even distance 
apart, which is of interest in its own right. 

A notable featnre of these processes is that they can be interpreted in 
a totally discrete manner. For finite graphs, they can be generated via a 
uniform probability measure on the possible orderings in which a finite col¬ 
lection of “events” associated with the vertices or edges of the graph can 
occur. For infinite graphs, they may be obtained by taking a limit of finite 
graphs. Moreover, blocking RSA and the annihilation process on graphs pro¬ 
vide a method of generating random independent sets of vertices, providing 
an alternative to the “hard-core” model which has received recent attention 
from discrete mathematicians [2, 9, 22]. 

2. Definitions and main results. 

2.1. Definitions for finite graphs. Let V be the vertex set, and E the 
edge set, of a graph (we consider only graphs with no loops or multiple 
edges). We start by defining onr models in the case where V is hnite. Each 
of them is described as a continuous-time stochastic process {XY,t > 0) 
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taking values in {0, l}'^, with XY {v) = 1 having the interpretation that site 
n G F is occupied at time t, and XY (v) = 0 having the interpretation that v 
is vacant at time t. For each v, the evolution of (XY{v),t > 0) is monotonic, 
so that the value of XY {v) changes at most once during the entire evolution 
of the process. 

Blocking RSA. Let the random variables {Ty,v £ V) be independent and 
uniformly distributed on (0,1), with Ty representing the arrival time of a 
particle at site v. All sites are initially vacant, and site v becomes occupied 
at time Ty if and only if none of its neighbors is already occupied. If site v 
has one or more occupied neighbors at time r„, the particle arriving at v is 
rejected (blocked) and there is no change to the state of the system at time 
Ty (this is the nearest-neighbor blocking/exclusion effect). Once occupied, a 
site remains permanently occupied. 

Dimer RSA. Suppose (re,e G E) are independent and uniform on (0,1). 
Here Te represents the arrival time of a dimer at edge e. All sites are initially 
vacant. 

Let e = {u,v} G FI be an edge of the graph, with endpoints u,v. If sites u 
and V are both vacant at time Te, then the dimer arriving at edge e at that 
time is accepted and sites u,v both become occupied at that instant. If sites 
u,v are not both vacant at time Te, then the dimer arriving at e is rejected, 
and nothing changes at time Te- Once occupied, a site remains permanently 
occupied. 

The annihilation process. Suppose (re,e G E) are independent and uni¬ 
form on (0,1). Here Tg represents the time of occurrence of an “attack” along 
the edge e. 

In this case each site is initially occupied. For each edge e = {u, u}, if sites 
u,v are both occupied immediately prior to time Te, then at time Te either 
site u or site v (each with probability 1/2) becomes vacant (i.e., one site 
“attacks” the other). If u,v are not both occupied at time Te, then there is 
no change in the state of the system at that time. Once it becomes vacant, 
a site remains vacant permanently. 

The multiple annihilation process (MAP). All sites are initially occu¬ 
pied, and {Ty,v G V) are independent and uniform on (0,1). At time Ty, site 
V (if still occupied) annihilates all of its neighbors and makes them all va¬ 
cant, while remaining occupied itself. If v is already vacant at time Ty there 
is no change to the state of the system at time Ty. 

Although defined in terms of annihilations, the MAP is most closely re¬ 
lated to the blocking RSA process. To demonstrate this relationship and 
aid understanding of the MAP, we introduce a further process {ZY ,t > 0) 
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with state-space {0,1, 2}^, which we call the undecided site proeess (or USP 
for short). Sites in states 0, 1, 2 are interpreted as vacant, occupied and 
undecided, respectively. In the USP, all sites are initially undecided (i.e., in 
state 2). If site v is undecided immediately prior to time then at that in¬ 
stant site V becomes occupied and all neighboring sites become (or remain) 
vacant. Otherwise, there is no change at time r^. 

With the undecided site process defined in this way, for u G U define 
coupled {0, l}^-valued processes {X^,t > 0) and {Y ^> 0) by 

(1) XY(v) = Y/iv) = ZY(v) if ZY{v) G {0,1}, 

( 2 ) XY(v) = 1 and yY{ v) = 0 ifZYiv) = 2. 

In other words, undecided sites count as occupied for the process {XY,t > 0) 
and as vacant for the process (YY ,t >0). It is not hard to see that with 
this construction, the process {XY ,t >0) is a realization of the MAP and 
(YY ,t >0) is a realization of the blocking RSA process. At the terminal 
time 1, there are no undecided sites left and therefore the ultimate state of 
each site is the same for the coupled MAP and blocking RSA process, that 
is, XY (v) = YY (v) for all v. In particular, the distribution at time 1 of the 
MAP and the blocking RSA process are identical. In view of this similarity, 
we shall treat the MAP only briefly (in Section 4.4). 

2.2. Remarks. For each of the blocking RSA, dimer RSA and annihi¬ 
lation processes, the evolution of {XY ,t > 0) is determined by the order 
of occurrence of the “event times” or Tg (together with the directions 
of “attacks” in the case of the annihilation process). Since they are as¬ 
sumed independent and identically distributed (i.i.d.), the “events” (arrivals 
of monomers, of dimers or attacks depending on the model) are equally 
likely to occur in any possible order. If we took the event times to be i.i.d. 
with some other nonuniform continuous distribution function F, then this 
would still be the case. In other words, changing the distribution function F 
amounts merely to a reparametrization of time. If the distribution F were 
chosen to be exponential, then the time-evolution of XY = {XY {v),v G V) 
would be that of a time-homogeneous continuous-time Markov chain with 
state-space {0,1}^. However, in fact it turns out to be more convenient to 
assume instead that F is the uniform distribution on the unit interval (0,1), 
in which case these processes run as non-time-homogeneous continuous-time 
Markov chains. 

Note also that while we could allow events to occur repeatedly at the 
same location (e.g., monomer arrivals at the times of a Poisson process for 
each site), in each model only the first “event” at any given location can 
possibly have any effect (e.g., if a monomer arriving at a particular site is 
blocked, then any later monomers arriving there are also blocked). Hence, 
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only the time of the first arrival is relevant, and as just mentioned, we can 
just as well (and do) take this to be uniformly rather than exponentially 
distributed. 

Owing to the monotonicity of the evolution, each process ultimately reaches 
a jammed state by time 1. For blocking (resp. dimer) RSA, jamming means 
that no vertex (edge) is any longer available to accept incoming particles. 

For annihilation, jamming means there is no pair of adjacent occupied sites. 
Both blocking RSA and annihilation produce a random independent set (in 
the graph-theoretic sense) of occupied sites at time 1; for blocking RSA 
the set of occupied sites is independent at all times. Of interest are the joint 
probability distributions of the states of sites at jamming, and also the states 
of sites at intermediate times prior to jamming. 

2.3. Extension to infinite graphs. Let us now consider our three main 
processes on a countably infinite graph. While it is not hard to envisage 
that the preceding definitions of our monomer, dimer and annihilation pro¬ 
cesses should carry through to the case where V is countably infinite, we 
need to check that these processes remain well defined. Evans ([7], page 
1302) appeals to general results from Liggett [11] to do this, but we prefer a 
more direct approach, since the results in [11] are stated for interacting par¬ 
ticle systems on rather than general graphs, and for time-homogeneous 
Markov processes rather than inhomogeneous ones, and in any case the di¬ 
rect approach is simpler here. 

A sufficient condition for our processes to be well defined on an infinite 
graph is obtained by assuming that V has bounded degree. In this case, the 
following result shows that we can define the blocking, dimer or annihilation 
process on E as a limit of processes on finite subgraphs. 

For any sequence of sets {Vn)n>i, we dehne the limit liminf(I4) := U^i r\m=n 

Proposition 1. Suppose {V,E) is a countably infinite graph with vertex 
set V and edge set E, with uniformly bounded vertex degrees. 

(i) Suppose Ti;,v €V, are independent and uniformly distributed on [0,1]. 
Then, except on an event of probability zero, it is the case that for all v G 

V and t G [0,1] and for any sequence {Vn)n>i of finite subsets of V with 
liminf(I4,) = P; for blocking RSA the limit 

(3) xY{v)-.= \\mjX^-{v)) 

exists and does not depend on the choice of {Vn)n>i- 

(ii) Suppose Te,e€ E, are independent and uniformly distributed on [0,1]. 
Then, except on an event of probability zero, it is the case that for all v G 

V and t G [0,1] and for any sequence {Vn)n>i of finite subsets of V with 
liminf(Ei) = V, for either dimer RSA or the annihilation process the limit 
(3) exists and does not depend on the choice of (Pn)n>i- 
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We use (3) to define the {0, l}^-valued process {Xt,t > 0), for countably 
infinite V. The almost sure convergence in (3) implies convergence of finite¬ 
dimensional distributions in the sense that that for any finite set of pairs 
, {tk,Vk) in [0,1] X V, and any sequence {Vn)n>i of finite subsets 
of V with liminf(W) = V, it is the case that 


{vfi, l<j<k)^ {Xl {vfi, l<j<k) 


as n —> oo. 


Both in the preceding result and throughout this paper, we are adopting 
the following notational conventions. Notation such as V,Vn ox W strictly 
speaking denotes the vertex set of a graph, but often will also be used to 
denote the graph itself. We often have occasion to consider one graph em¬ 
bedded in another; when V <ZW and W is the vertex set of a graph, then 
(e.g., for RSA on V) we assume the graph structure on V is that induced by 
the graph structure on W. In other words, two elements of V are adjacent 
(in the graph V) if and only if they are adjacent in the graph W. 

Another way to define the process (X^,t > 0) on infinite graphs is by a 
recursion which we illustrate only for the case of blocking RSA. If (Ty,v & V) 
has a local minimum at site w, then denote site w an aceepted site and denote 
all neighbors of w as rejected sites. Remove all accepted and rejected sites 
from V and repeat. Continue indefinitely to generate a set of accepted sites. 
For each site v, and t > 0, set XY (v) = 1 if and only if t > and site v is 
accepted; otherwise set XY (v) = 0. It can be shown that if V is countable 
with bounded degree, the process XY generated this way is identical to 
that given in Proposition 1. For a continuum blocking RSA model, a similar 
approach is used and described in more detail by Penrose ([18], page 160). 

2.4. Forward equations. Given a vertex set P of a graph, write dist(u, w) 
for the graph distance between two elements of V. We introduce the notation 


Afk{v, V) := {w : dist(u, w) < A:}, k> 0 


( 4 ) 


for neighborhoods of vertices in V, and write just Nkiy) for J\fk{v,V) if the 
graph V is clear in the context. 

The key to quantitative results is the fact that we can obtain differential 
equations for the probability of vacancy at particular sites. We hrst give 
such a result for blocking RSA. 

Proposition 2. Let V be finite or countably infinite with bounded de¬ 
gree. Let W be a finite subset of V. Then for 0 <t <1, for blocking RSA, 


( 5 ) 


±P[xY{v) = 0,v€W] 

= -J2 P[Wi'"''^"^(u;)=0,u;G(WUA4(i^,P))\W]. 


v&W 
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We turn to the corresponding differential equations for dimer RSA and 
the annihilation process. Consider first the special case of vertices of degree 

1 . 


Proposition 3. Let V be the finite or countably infinite vertex set of a 
connected graph with bounded degree and more than one edge. Let W be a 
finite subset of V, and suppose each element of W has degree 1. Then for 
0 < t < 1, for dimer RSA, 

±p[xY{v) = o,veW] 

( 6 ) 

= - X! = o,iuG (wuM(i^,'C))\ W] 

vew 

and for the annihilation process, 

±p[xY{v) = i,veW] 

(7) 

= - i E = l,we{WuAf,{v,V))\ {u}]. 

^ v£W 

The next result shows that the general case can be reduced to the case 
just considered, of vertices of degree 1, by considering a modified graph Gw, 
which is illustrated in Figure 1. 

Given a graph G = {V,E), and given a vertex w oi degree d{w), 
enumerate the vertices in Mi{w, V) as ni{w),... and let = 

be the graph that is obtained from G by replacing the vertex 
w with d{w) vertices of degree 1, denoted ujf,..., with v* connected 

to nfiw) by an edge for 1 < z < d{w). 

Given a finite subset W of V, let Gw = {Vw,Ew) be the graph obtained 
from V by successively replacing each vertex tc in VF by a collection of d{w) 
vertices of degree 1 as described above (the order in which the vertices in W 
are taken is immaterial). In other words, define Gw for W CV of cardinality 
k, inductively for A: = 1,2,3,..., by 

( 8 ) Gwu{w}iP) = {Gw){w}- 

For finite W C I/, let W* be the set of vertices in Vw winch replace the 
original vertices in W, a total oiJ2wewd{w) vertices. Both Gw and W* are 
illustrated in Figure 1. 

Proposition 4. Let G = {V, E) be finite or countably infinite with bounded 
degree. Let W be a finite subset of V. Then for 0 <t < 1, for dimer RSA, 

(9) P[XY{v) = 0,veW] = PIXY'^ {v) = 0, u € W*] 
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Fig. 1. An example of the construction of Gw ■ Here W is the set of white vertices in 
the graph G on the left, and W* is the set of white vertices in the graph Gw on the right. 


and for the annihilation process, 

(10) P[XY{v) = l,veW]= (-y) = l,veW*]. 

2.5. Summary of exact results for regular trees. Stochastic processes on 
trees have been the subject of much activity in the probability community 
in the last 10-20 years. See, for example, [12, 17]. One reason for looking at 
such processes is that they are sometimes more tractable than the equivalent 
processes on lattices such as while providing some insight into the latter 
processes. 

In Section 4 we consider our RSA and annihilation processes on the {k + 
1)-regular tree (Bethe lattice), which we denote by T^. This is a countably 
infinite connected graph in which each vertex has degree A: -|- 1 (A: > 1 a hxed 
integer). We shall give exact computations of probabilities on Tk which we 
summarize as follows. 

For an arbitrary vertex u of T^, let a{t) denote the probability that vertex 
V is still in its initial state (vacant in the case of blocking or dimer RSA on 
Tk, occupied in the case of annihilation on Tk) at time t. Let /3(A) be the 
probability that v is still in its initial state at time t for these processes on 
a graph obtained from Tk by removing one branch from v (in the case of 
blocking RSA) or by removing every branch except one from v (in the case of 
dimer RSA or the annihilation process). We shall use the forward equations 
from the preceding section to derive the following formulae for k>2 (only 
the second line needs modification for A: = 1): 


Blocking RSA Dimer RSA Annihilation 


P{t)=-f){tf 

a'{t) = -P{tf+^ 
a{t)^\{l + (5{tf) 


f}'{t) = -fj{tf 

a{t)=m'‘^^ 


/3'(t)=-/l(t)V2 


rp 

yielding the following numerical values for the occupation probability P[X^ (u) = 1] 
at time A = 1: 
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Occupation probabilities at time t = 1 on Tfe 

k 

Blocking 

Dimer 

Annihilation 

1 

0.432 

0.865 

0.368 

2 

0.375 

0.875 

0.296 

3 

0.333 

0.889 

0.250 

4 

0.302 

0.901 

0.217 

5 

0.276 

0.911 

0.192 

10 

0.200 

0.940 

0.124 

20 

0.135 

0.964 

0.074 


The formulae above for a{t) have previously been seen in the literature, 
not always in this form. In early work, Widom [23] gave an exact formula 
for the occupation density for blocking RSA in T 2 but later retracted this 
[24], The correct formulae are given, at least for t = 1, by Evans ([5], page 
2529) for dimer RSA, by Fan and Percus ([8], equation (4.8), see also [6]) for 
blocking RSA and by O’Hely and Sudbury [14] for the annihilation process. 
However, the first two of these derivations are quite involved, and proba- 
bilists may find our derivation more appealing. 

We shall also give exact formulae for the covariances between the states 
of different sites of T^, for each model. The formulae are given in Section 4, 
and in this summary we just give some of the resulting numerical values of 
correlations between two sites at time t = 1. Clearly these depend only on 
the graph distance between the two sites (and on k): 


Correlations 

Distance 

Blocking 

Dimer 

Annihilation 

fe = 3 

fe = 5 

fe = 3 

fe = 5 

fe = 3 

fe = 5 

1 

-0.5000 

-0.3820 

-0.1250 

-0.0982 

-0.3333 

-0.2383 

2 

0.1760 

0.1003 

0.0070 

0.0103 

0.0758 

0.0407 

3 

-0.0474 

-0.0200 

0.0085 

0.0007 

-0.0130 

-0.0054 

4 

0.0103 

0.0032 

-0.0041 

-0.0005 

0.0018 

0.0006 

5 

-0.0019 

-0.0004 

0.0012 

0.0001 

-0.0002 

-0.0001 


See [16] for previous work on time-evolution of correlations for RSA on Ti. 
We believe that our formulae for correlations in Tfc,fc > 2, are new. 

We note that the blocking and annihilation processes, but not the dimer 
process, show a pattern with sites at odd distances being negatively cor¬ 
related and those at even distances being positively correlated. Next we 
note that correlations fall off extremely fast, since they are truncations of 













10 


M. D. PENROSE AND A. SUDBURY 


the exponential series (see Section 4). Also, correlations for the annihilation 
process are less than those for the blocking process. This is becanse for one 
site to influence another, not only do the “events” along a path between 
the two sites have to be in the correct temporal order, in the case of the 
annihilation process they have to be in the correct direction as well. Finally, 
the correlations for k = b are smaller than those for k = 3. 

2.6. Bounds for bipartite graphs. A graph is bipartite if its vertex set can 
be partitioned into two independent (in the graph-theoretic sense) subsets. 
The next result says that for either RSA or the annihilation process, on a 
bipartite graph, the states of sites an even graph distance apart are positively 
correlated and the states of sites an odd graph distance apart are negatively 
correlated. This shows, among other things, that the pattern of positive and 
negative covariances for blocking RSA or the annihilation process, seen in 
the table from Section 2.5, extends to all k and all distances (since any tree 
is bipartite). 

Theorem 5. Suppose V is the vertex set of a bipartite graph. For j = 
0,1,..., fc, suppose Wj C V, with Vj € Wj and tj € [0,1], and define the event 
Ej by 

Ej := (vj) = 0} if dist(nj, no) is even, 

Ej := (vj) = 1} if dist(nj, no) is odd. 

Then for either bloeking RSA or the annihilation proeess, 


(11) 

p 

k 

nu 

> P[Eo]P 

k 

flu 



Lj=o J 


-i=i 


In Theorem 5, it is to be understood that the processes {X^ft > 0) are 
defined for all j in terms of the same set of random arrival times (r^,n G 
V) (in the case of RSA) and the same set of random arrival times Tg and 
random directions of attacks, along each edge (in the case of the annihilation 
process). The main case of interest in Theorem 5 is when Wj = V for each 
j, but we shall also use the result in the more general form given. 

Let us say R is a bipartite lattice if (i) V is the vertex set of a bipartite 
graph and (ii) for any vi,V 2 ,wi,W 2 €V with vi adjacent to V 2 and rci adja¬ 
cent to W 2 , there exists a graph isomorphism from V to itself which sends Vi 
to Wi, i = 1,2. Examples include the integer lattice (with edges between 
each pair of vertices a unit Euclidean distance apart), and also the hexag¬ 
onal lattice (but not the triangular lattice which is not bipartite). Various 
finite graphs are also bipartite lattices, including any lattice torus and also 
the discrete n-cube. 






DEPOSITION PROCESSES ON GRAPHS 


II 


For a bipartite lattice, every vertex has the same degree. One of our main 
results says that for blocking RSA on such a lattice, the exact values given 
in the table for vacancy probabilities in the corresponding tree with the 
same vertex degrees, provide upper bounds for vacancy probabilities for the 
bipartite lattice. 

Theorem 6. For blocking RSA on a bipartite lattice with vertices of 
degree k + 1 where k>2, for any vertex v we have 

<-(l+ (i^(^_l)^)2/(fc-i))- 

In partieular, 

P[XY(v) = 0] < i(l + so P[XY(v) = 1] > i(l - 

In the special case of the square lattice we have k = 3 so that by 
Theorem 6, P[XY(v) = 0] < the occupation probability at jamming 

satisfies P[XY(v) = 1] > 1/3. Simulations suggest a value of approximately 
0.364 ([7], page 1292). 

In the case of the hexagonal lattice, we have k = 2 so that by Theorem 
6, P[XY(v) = 1] > 3/8, which compares favorably to the lower bound of 
0.3596 provided (with a good deal of effort) by Caser and Hilhorst in [3]. 
Simulations indicate a value of 0.379 [6, 7, 23]. 

For the integer lattices and we have k = 5 and k = 7, respectively, 
so that Theorem 6 yields lower bounds for P[XY (v) = 1] of 0.2764 and 
0.2386, respectively, compared to simulated values [7, 13] of 0.304 and 0.264, 
respectively. 

Remark. Our attempts to obtain good rigorous bounds for the annihi¬ 
lation process on bipartite lattices, using a similar argument to the proof of 
Theorem 6, have not been so successful. See the remarks at the end of the 
proof of Theorem 6. 

2.7. Results on normal approximation. The next result is a normal ap¬ 
proximation result for the total number of occupied sites at time t on some 
(large) finite set V for one of our processes on a (finite or infinite) graph 
with vertex set W satisfying V C W. We denote this quantity St{V, W), that 
is, set 

(12) St{V,W):=J2xY^{v). 

vev 

Let <h(-) denote the standard normal distribution function, and let card(F) 
denote the number of elements (cardinality) of a finite set V. 
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Theorem 7. Let D >2 be a finite constant and let e > 0. Then there 
exists a finite constant C = C{D,e) such that for any finite V, any finite 
or countable W such that V C W, and any graph with vertex set W and 
maximum degree at most D, and all x G M, we have 


(13) 


f St{V,W)-ESt{V,W) \ ^ ^ 


$(x) 


< C Var (1/, IT)) ■ (card(l/)) 


for any of blocking RSA, dimer RSA or the annihilation process. 


Suppose we are given a sequence of pairs (14,, where for each 

n, Wn is the vertex set of a finite or countable graph and 14 4 Wn, and 
card(14) —> oo as n —> oo. For example, 114 could be the same as 14 for each 
n, or all 114 could be the same infinite vertex set IF, containing all the 14- 
Under these circumstances, Theorem 7 gives us a central limit theorem with 
error bound on the rate of convergence, provided the Wn have uniformly 
bounded degrees and we have reasonably fast growth of Var(S't(14, Wn)). 

It can be shown that Var:{St{Vn,Wn)) = 0(card(Fi))j and one might ex¬ 
pect Var(S't(14, Wn)) = 0(card(Fn)) which would be sufficient to guarantee 
a central limit theorem from (13) with a rate of convergence card(14)^~^'^^- 
The next resnlt guarantees that this is indeed the case for t < 1. 


Theorem 8. Let t G (0,1). Let D < oo. Then for any F C IF with all 
degrees in W bounded by D, for blocking RSA, 

(14) Var(5t(F,lF)) >t(l-t)^+H^ + l)”^card(F). 

If there are no isolated vertices in V, then for either dimer RSA or the 
annihilation process, 

(15) Var(St(F, IF)) > (l/2)t(l - t)'^^-^{2D - 1)"^ card(F). 


In the case t = 1, certain counterexamples suggest that, at least for block¬ 
ing RSA or the annihilation process, there is no very simple condition to 
ensure that Var(S'i(F, IF))/card(F) is bounded away from zero. These are 
illustrated by Figure 2 and described below. 

For the first example, suppose Vq is an arbitrary finite graph, and V is 
obtained by starting with Vq and adding a “twin” vertex for each vertex in 


I I I 1 I 1 I 1 I vvvvvvvvv 

Fig. 2. Two graphs for which the ultimate number of particles deposited has zero vari¬ 


ance. 
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Vo, with each added vertex connected by an edge only to its twin. Ultimately, 
for each pair of twinned sites exactly one of the twins will be occupied, so that 
the variance of S'i(U, U) is zero, for either blocking RSA or the annihilation 
process, although V can be arbitrarily large. 

For a second example, suppose Vq is an arbitrary finite graph, and V is 
obtained by starting with Vq and adding for each vertex v G Vq a pair of 
extra vertices v',v", with u' connected by edges only to v'' and to v, and v" 
connected only to v and to u'. In other words, a triangle is appended to each 
vertex in the original graph Vq, and again we conclude that the variance of 
5'i(U, U) is zero for either blocking RSA or the annihilation process. 

The first counterexample suggests that we may be in difficulties trying 
to find lower bounds for Var(5'i(U, U)) if all of V lies near to the boundary 
(i.e., to the set of vertices of degree 1) while the second suggests that there 
can be problems if V is not bipartite. 

To give a sufficient condition for (14) to hold even for t = l, we introduce 
further terminology. Let us say that a site u in a graph W has positive entropy 
if (i) there exist at least two sites in W that are adjacent to v but not to each 
other and (ii) there exists an independent subset of W) \Af 2 {v, W) 
such that every vertex in M 2 {v, W) has at least one neighbor in B^. We shall 
later call B^ a blocking set for v. 

A sufficient condition for v to have positive entropy is that W be the 
vertex set of a bipartite graph and u be a graph distance at least 3 from the 
boundary of W [one can then take B^ = Mz{v^W)\M 2 {v.,W)]. 

Denote the set of sites in W with positive entropy by W~^. 

Theorem 9. Let D G N. For any graph W and any V CW with all 
degrees in W bounded by D, for blocking RSA, 

(16) Var(5i(U,lU))> (^^^^^^^^^)card(lU+nU), 

and likewise there is a constant r] > 0 depending only on D such that for the 
annihilation process, 

(17) Var(5i(U,iy)) >r 7 card(W+nU). 

It should be possible to find an analogous result to Theorem 9 for dimer 
RSA, but we have not done this. 

We can easily deduce the following central limit theorem from the pre¬ 
ceding results. Let A7(0, cr^) denote the normal distribution with mean zero 
and variance 

Theorem 10. Suppose {Wn)n>i is a sequence of finite graphs with uni¬ 
formly bounded degrees, and 14 , is a sequence of subgraphs of 114 and with 
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lim„_,oo card(V^) = oo. Then for 0 < t < 1, 


(18) 


StiVn^Wn) - ESt{Vn,Wn) V 


AA(0,1) 


^\aTSt{Vn,Wn) 

for any of blocking RSA, dimer RSA or the annihilation process. If, in ad¬ 
dition, the sequence of pairs {Vn,Wn) satisfies 

lim inf (card (14 H card(14)) > 0 , 


then (18) also holds for t = 1 in the case of either blocking RSA or the 
annihilation process. 


Theorem 10 has some overlap with Theorem 2.2 of [19]. However, Theorem 
10 addresses processes on general graphs whereas [19] restricts attention to 
certain graphs embedded in Z'^. Also, the proof of Theorem 10 provides a 
bound on the rate of convergence, unlike the corresponding proof in [19]. On 
the other hand, [19] is concerned with a more general class of processes than 
those considered here, and provides information on convergence of rescaled 
variances (for the graphs covered by [19]). 

Finally, we can give more explicit central limit theorems in the one¬ 
dimensional case. Let Z be endowed with the usual nearest-neighbor graph 
structure (so that it is isomorphic to Ti from Section 2.5). Let denote 
the set {1,2,..., n}, that is, a line of n adjacent sites in Z. 

Theorem 11. Let t G (0,1], and let (y‘^\^{t) :=te~^^. For blocking RSA, 


as n —> 00 


(19) 

n ^Yav St{Ln,T,) ^ a^i^{t), 

(20) 

n~^YaTSt{Ln,Ln) 

(21) 

n-\St{Ln,Z) - ESt{Ln,Z)) ^ N{t),al^{t)) 

and 


(22) 

n-\St{Ln,Ln) - ESt{Ln,Ln)) 


Moreover, for dimer RSA and for the annihilation process the results (19)- 
(22) still hold with the limiting variance (T3L replaced by cr^j and re¬ 
spectively, given by 

'^Di(^) : '^an(^) (1 T 2t)e — e 

We note that using different methods, (19) was previously obtained, for 
blocking RSA only, by Pedersen and Hemmer [16], and the special case t = l 
of (20), for dimer RSA only, was obtained long ago by Page [15]. 
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3. Proof of existence results and forward equations. 

Terminology. By a path of vertices in a graph we mean a sequence of 
distinct vertices such that any two consecutive vertices in the sequence are 
connected by an edge. By a path of edges, we mean a sequence of distinct 
edges such that any two consecutive edges in the sequence have an endpoint 
in common. 

Proof of Proposition 1. (i) Suppose W dW ClV, with W and W 

both finite. Then for v G W, it is the case that {v) = (v) unless 

there is a path of vertices vo,vi,... ,Vm, starting at vg = v and ending at 
some vertex Vm G W \ W, such that > Ty^ > ■ ■ ■ > since the only 
way that the effect of a vertex in W' \ W can propagate to v is along such 
a path in reverse order. 

Therefore, for infinite V, a sufficient condition for the existence of the 
limit (3) for u G P is the existence of a finite random variable R{v) such 
that there is no path {vg, vi,..., Vm) in V with vg = v, m> R{v), and with 
Tvq> Tvi > ■ ■ ■ > Tv^, since in that case, the value of X^ (v) will be the same 
for all finite W <ZV containing all vertices within graph distance R{v) of v. 

Let V gV. Since degrees are assumed bounded by a constant which we 
denote D, the number of distinct self-avoiding paths starting at v of length 
m is bounded by D”*. For any such undirected path (uq, vi,..., Vm), we have 

P[tvo >Ty,> -> = l/(m + 1)! 

since any ordering of the times Ty ., 0 < j < m, is equally likely. Hence, the 
probability that there is a path of length m starting from v with associated 
arrival times in decreasing order is bounded by D'^/{m + 1)!, and by the 
Borel-Cantelli lemma, the total number of paths starting from v with asso¬ 
ciated arrival times in decreasing order is almost surely finite; we can take 
R{v) to be the maximum of the lengths of such paths. This goes for any 
V gV , and the existence of the limit in (3) for all v,t follows. 

(ii) Given W dW dV, with W and W hnite, and given v G W, it is the 
case that the value of X^ (v) = X^ (v) for dimer RSA, unless there exists 
a path of edges ei,e 2 ,...,em that starts at v and ends at some vertex in 
W' \ W, such that Tei > Teg > ■ • • > Te^ ■ 

We assert that there exists almost surely a random variable R'(v} such 
that there is no path of length greater than R'{v), starting at v, with all 
edge arrival times along the path in decreasing order. The proof of this is 
the same as for the corresponding assertion in the proof of part (i). 

Then X^ (v) will be the same for all finite W dV with W containing 
all vertices a graph distance at most R'{v) from v. Hence, (3) holds for the 
dimer process. 
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Next consider the annihilation process. With v, W, W' as before, this time 
we have (v) = (n) unless there exists a path of edges ei, 62 ,..., Cm 

that starts at v and ends at some vertex in W' \ W, such that Te^ > Tej > 

• • • > and in addition all “attacks” along the edges are in the direction of 
the path towards v. Hence with R\v) as defined for dimers, the annihilation 
process X^ (v) will be the same for all finite W with W containing 
all vertices a graph distance at most R'{v) from v. Hence, (3) holds for the 
annihilation process. □ 


Proof of Proposition 2. Suppose V is hnite or countably infinite 
with bounded degree, and VP is a finite subset of V. For v and 0 < t < 
t + h<\, we have 

P[{XY {w) = 0,n; € W} n = 1}] 

= P[{tv G (t,f + h]} n {xY (w) = 0,w e {W Uj\fi{v)) \ {n}}] + o{h) 

= hP[xY {w) =0,w e {W U A/i(n)) \ {n}|rt, G (t,t + h]\ + o{h) 

= hPlxY^^""^ {w) = 0,ui G (VF U \ M] + o{h). 

Summing over v €V we obtain 


(23) 


{xY{w)=0,w€W}nl \J{Xt+Uv) = '^} 

= hlj2 ^’[^r^^'’^(w') = 0,u;G(VFuA4(ii))\{i^}]) +o{h). 
\t;eVF / 


Similarly, for 0 < h <t, 


{xY_^,iw)=o,wGW}n( U{W(i;) = i} 


^v€W 


(24) =h('^P[X,':\^">(n;) = 0,n;G(VFuM(i;))\W])+o(/i) 

\t;eVF / 


= /i( ^ P[xY^^^\w) = 0,w e (VF U A/i(ii)) \ {'!’}] I + o{h), 

Vi;GVF / 

where the last line follows from left-continuity in t of P[xY^^^\w) = 0,ii; G 
(VFUA/i(ii))\{t’}], which is easily seen. The result (5) follows from (23) and (24). 

□ 


Proof of Proposition 3. Let V be the finite or countably infinite 
vertex set of a connected graph with bounded degree and more than one 






DEPOSITION PROCESSES ON GRAPHS 


17 


edge. Let W he a finite subset of V, and suppose each element of W has 
degree 1. Let v € W. Then Mi{v) consists of a single vertex, denoted u, 
which is not itself in W (otherwise, G would either be disconnected or have 
only a single edge). Then for dimer RSA, writing e for the edge {u,v}, we 
have 

P[{XY (w) =0,w GW} n (v) = 1}] 

(25) = P[{Te G{t,t + h]} n {X^(w) =h,w g{W U A/i(u)) \ {u}}] + o{h) 

= hP[X^^^'’^ {w) = h,w G{W uNi{v))\{v}]+ o{h). 

Summing over v gW , we obtain 


P 


(26) 


{xY{w) = h,wGW}r\( U {^i+/.(^) = 1} 




and by arguing as for (24) we also have a similar expression to (26) for 


= hY, P[X^^^''\w) = 0,wG{WuMi{v))\{v}]+o{h), 

v£W 


P 


{xtH{^) = o,w<^w}n 


U = 1} 

vew 


Thus, we obtain (6) for dimer RSA. 

For the annihilation process, given the event of an attack occurring along 
edge {u,v} between times t and t + h, this attack has conditional probability 
1/2 of being in the direction from u to v, so that instead of (25) we obtain 

P[{XY (w) = l,wGW}n {xY+dv) = 0}] 

= {h/2)P[xY'^^'’Hw) = hwG{WUAfi{v))\{v}]+o{h). 

The rest of the derivation of (7) for the annihilation process is the same as 
the derivation of (6) for dimer RSA. □ 


Proof of Proposition 4. Let IT be a finite subset of V. Th^e is an 
obvious canonical bijection between the edges of G and the edges of Gw (see 
Figure 1), and using this bijection, for either the dimer RSA or annihilation 
process there is a coupling between a realization of the process {XY ,t > 0) 

and the process {xY '^, t > 0) where we take the arrival time to be the 
same for edge e of G and for the corresponding edge of Gw (and in the 
annihilation model, we take the direction of attacks to be the same in both 
processes). 

With this coupling, the state of each vertex v gV for the process XY is 
the same as the state of the corresponding vertex or vertices in Vw for the 
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process , up to the instant when some site v first changes state, at 
which point one may see a difference between the coupled processes because 
only some of the sites in VV change their state. 

Thus, with this coupling of processes, the first instant at which some site 
in W changes its state in the process ) is the same as the first instant at 

which some site in W* changes its state in the process (X^^), and therefore, 
regardless of the coupling, the first instant at which some site in W changes 
its state in the process (X^) has the same distribution as the first instant 

at which some site in W* changes its state in the process {X^^)- This is 
true for both dimer RSA and the annihilation process, so that (9) and (10) 
hold. □ 


4. Analysis of the fc-tree. This section is concerned with deriving the 
results for the case where V is the vertex set of the (k + l)-regular tree T^, 
which were summarized in Section 2.5. 

It is clear that for each of the processes under consideration, for vertices 
v,w of Tfc the covariance of xj*‘{v) and xj’^{w) depends on v,w only via 
dist(n,r(;). We shall use the notation 

(27) Cov{m,t,Tk)'■=Coy{xJ*‘{v),xJ^{w)) if dist(n,r(;) = m. 

Naturally Cov{m,t,Tk) will also depend on whether we are considering the 
blocking RSA, dimer, annihilation or multiple annihilation process. 

Let be the rooted {k + l)-tree, that is, the graph with a single vertex 
(the root, denoted vq) of degree k and all others of degree k + 1. 

4.1. Blocking RSA on T^. 


Theorem 12. For blocking RSA on T^, it is the case that for arbitrary 

v^Tk, 

(28) P[Xf^{v) = 0] = i(l + {P[xf^{vo) = 0])2), 


with 


(29) P[xf^ivo) = 0] 
In particular, 




k>2, 
k = l. 


P[xf^{v) = 0] 


|(l + /c-2/(^-b), k>2, 


(30) 


k = l. 
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Proof. Fix k, and define f3{t) and a{t) (for an arbitrary vertex v G Tfc) 
by 

(31) P{t) = P[xf^vo) = 0], a{t) :=P[X^^{v)=0] 


(clearly the second probability does not depend on the choice of v). The 
subgraph of induced by removing uq consists of k unconnected copies of 
T^, so by Proposition 2, 

( 32 ) p'it) = 


with initial condition /3(0) = 1, together giving us /3(t) = e * for k = 1, and 


|(/3(t)'-")=A:-l ^ /3(i)i-fc = i + (A:-l)t, 

and (29) follows. 

In a similar manner to (32), we obtain 


(33) 


a'{t) = 




d 

dt\ 2 J 


k>2, 


with q:( 0) = /3(0) = 1, so that a{t) — 1 = (/3(t)^ — l)/2 and hence a{t) = 
(/3(t)^ + l)/2, that is, (28) holds. □ 


We now turn our attention to correlations between the states of two sites 
in Tfc. Recall the definition (27). 


Theorem 13. Let k >1 be an integer, and consider blocking RSA on 
Tfc. Set (5{t) := P[x'^'^ (uq) = 0]. Then for 0 < t < 1, and for integer m>l, 


(34) 


Cov{m,t,Tk) = 


(3{tf ^ (21og/3(t))"^+i+2» 
2 (m + l + 2n)! 


Remark. In the case k = \, using (29) we see that this agrees with the 
formula of Pedersen and Hemmer [16] (our time parameter t corresponds to 
1 — in their parametrization of time; see our remarks in Section 2.2). 


Remark. In the case where dist(u,r(;) = 1, our formula for the covari¬ 
ance in Theorem 13 simplifies to —(1 — a(t))^, with a{t) given in (31). This 
can be seen directly since for monomer RSA on any graph, P[XY (u) = 
XY (w) = 1] = 0 for adjacent sites v,w. 
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Proof of Theorem 13. Again define a{t) and /9(t) by (31). Set 

u{t) = 2\og{f3{t)). 

If k>2, let T” denote an infinite tree in which a single vertex (the root, 
denoted vq) has degree k — 1 and all the other vertices have degree k + 1. If 
k = l, let T^* consist of a single isolated vertex. Set 

7(t) :=P[xf“(^o)=0]. 

Then by Proposition 2 and (32), 

(35) 7'(t) = -/3(i)"-' = «'(i)/2. 

Since 7 ( 0 ) = 1 and u{0) = 0, we obtain 

(36) ^{t) = l + u{t)/2. 


For integer m > 1, let denote an infinite tree in which two vertices 
(denoted fo,!*!) have degree k and all other vertices have degree k + 1, and 
in which dist(no,r’i) = m. Set 

am{t) ■= (v) = 0 and Xj'’" (w) = 0], v,w€Ti^, dist{v,w)=m, 

Pm{t) '■= P\x'^*" (t^o) = 0 and xj'" (w) = 0], re G T^, dist(?;o, ■;«) = m, 

7m(i) := P[xf^’^{vo) = 0 and = 0]. 


Also set ao(t) = a{t), Po{t) = (3{t), 7 o(t) = 7 (t)- Then for m > 1, Proposition 
2 followed by (35) yields 




-2(3{tf Sm-l(i) 


d'lm 

du 


lm-li 


with 7 m( 0 ) = 1 and u(0) = 0, and with 7 o(t) = 1 + u{t)/2 by (36). This 
system of equations has solution 


(37) 


1 ™ 


7m = 


2 (m + 1 )! 


+ E 

Tl=0 


ni 


I ’ 


so that 7 m(t) ^ as m —> 00 , as would be expected. 

Now consider Pmit). By Proposition 2 and (35), we have for m > 1 that 

with /3m(0) = 1 and u(0) = 0. This system of equations has solution 


1 / m ^ 

1 P U 


/5m - X 1 + E 




n=0 


ni 
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and as m —> oo we have jdm —> |(1 + = a{t)P{t), as expected. 

Similarly for am{t), by Proposition 2 and (35) we have 

(38) aUt) = -2P{t)’^Pn,-i{t) ^ = 

au 

with am(0) = 1 and u{0) = 0, giving us 


l + (-l)™ e“ 
; I TT 


/ L(m-1)/2J ^m-l-2n ' 


Hence, 


1 e“ 

lim a 2 m = 71 + — (1 +sinhn) 

m^oc 2 2 


+ 2e“ +1 I \ T 

=-^-= —(l + coshn)= hm a 2 m+i 

4 2 m—»oo 

SO that limm^oo Oim{t) = as one would expect. This implies the covari¬ 

ance between sites a distance m apart is 


2 oo 

Cov{m,t,Tk) =am{t) - a{tf = -xe“ ^ 


,m-\-l-\-2n 


2' ^^{m + l + 2n)\' 


which implies the asserted formula (34). □ 


4.2. Dimer RSA on T^. In this section we derive the results for dimer 
RSA analogous to Theorems 12 and 13 from the preceding section, namely, 
(39) and (45) below for k>2, and (40) and (46) for k = l. 

For dimer RSA on T^, and for an arbitrary vertex u of T^, set 

a{t) = P[X^^ {v) = 0] , Pit) = Pixf^ (no) = 0], 

where has vertex vq of degree 1, all other vertices of degree A: + 1. Then 
the graph (Tfc){.y} described at ( 8 ) has A: + 1 components, each of them 
isomorphic to T^. By Propositions 4 and 3, 

ait) = Pit)>^+\ (3'it) = -Pit)\ 

and the solution to these [with P{0) = 1] is 
/3(i) = (l + (A:-l)t)-V(fc-i), 

(39) 

a{t) = (1 + (A: - i)t)-(^+i)/(fc-i)^ > 2 , 

P{t) = e~^, a{t) = e~'^^, 


(40) 


A: = l. 
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Set 

«m(i) := {v) = 0 and XP {w) = 0 ], v,w ^ T^, dist(n, w) = m. 

Define to be the tree with vertices no,ni of degree 1, with dist(no,ni) = 
m, and all other vertices of degree k + 1. By Propositions 3 and 4 we have 
am{t) = where 


/3^(t) = P[xf^-(no) = = 0], 

Then (3i{t) = 1 — t, and by Propositions 4 and 3, 

(41) /3:,(t) = -2/3(t)"-Vm-i(t). 

To solve these, set v{t) := 21og/9(t). Then v'{t) = —2(5{t)^~^ and 
( df3m Pmi't) Q 

with /3m( 0 ) = 1 and with 

/c 

/3l = 1 — t = - - -- - -; 


Ti 


m>2, 


k>2. 


k — 1 k — 1 

For any constants a, b, c, with c 7 ^ 0, the solution to the system of equations 
(42) with /3i(n) = a — is 


(43) 


Pm = 


YY) 1 TTX —2 n 

av 


E u 

71“ 


— j! 


b^m-l g (v/cV 




j=m-l 


In the present case (for k > 2) we have a = k/{k — 1), b = l/{k — 1) and 
c = 2/{l — k) so that 


m—1 


Pm{t) = 


fc(21og/3(t)) 

(m — 1 )!(A: — 1 ) 


^ /^^_(21og/3(t))^^i 

V i =0 


3 


i! 


) 


+ 


1 


2 Vl - A; 


- “ {{\-k)\ogP{t)p 

7! 

j=m—l * 


and hence for k>2, 


Cov{m,t,Tk) 

= am{t) - a{tf = P{tf^{Pm{t) - P{tf) 


= Pit) 


2 k f A:(21og/3(t)) 


m—1 


(45) 


y (m — l)!(fc — 1) 

/(21og/3(t))^' 1 


- E 

j=m-l 




2\l-k 


iil-k) logPit))^\ 
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For /c = 1 we still have (41) with /3i(t) = 1 — t, but now we have P{t) = e“*, 
so that f3{t) (for dimer RSA on Ti) equals 7m-i(t) (for monomer RSA on 
Ti, as given in Section 4.1), because they satisfy the same equations and 
agree for m = 1. Thus 


Pm{t) 


1 {-2t)^ ^ (-2t)’" 

2 ml “1 n! 


and hence, since (44) still holds for k = l, 


(46) 


Cov(m, t, Ti) 


—e 


-2t 


1 {-2t)^ 

2 ml 


OO 

+ E 

rx=m+l 


(-2t)M 

n! J ■ 


By (45) and (46), for all k we have Cov(l, l,Tfc) = —a(l)^ which can also be 
seen directly since for adjacent sites u, v we have P[XY (u) = XY (v) = 0] = 0 
for dimer RSA. 


4.3. The annihilation process on Tk- In this section we derive results for 
the annihilation process analogous to Theorems 12 and 13, namely, (47) and 
(49) below for k>2, (48) and (50) for k = 1. For the annihilation process, 
with as defined in the preceding section, let 

a{t) = P[xY>‘iv) = 1 ], Pit) = P[xf^ivo) = 1 ]. 

Then Proposition 4 gives us a{t) = P{t)^~^^, and Propositions 3 and 4 give 
us P'it) = —P{t)^/2. Also, /3(0) = 1. Hence, 

Pit) = il + ik-l)t/2)-^/^^-^\ 

(47) 

ait) = (1 + (A: - l)t/2)"^^+^)/(''“^\ k>2, 

(48) Pit)=e~^^‘^, a(t) = e“*, A: = l. 

Let T^^ be as described in the preceding section, and let 

o-mit) ■= P[xY'’ iv) = Xj^ iw) = 1 ], dist(u, w) = m, 

P^it) := P[xf^’-ivo) = xf^’-ivY = 1 ]. 

Then amit) = Pit)"^^Pmit), Piit) = 1 — t, and by Propositions 3 and 4, 

= Pm-lit), m>2. 

Let vit) = 2\ogPit). Then v'it) = —/3(t)^“^. Hence, 

dPm 


dv 


Pm—l, 
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with 




2g(l-fc)i;/2 

k-1 


k>2. 


Using the general solution (43) with a = {k + l)/{k — 1) and 6 = 2/(fc — 1), 
and c = 2/(1 — k), for k>2 we obtain 


^ (fc + l)(21og/l(t))—1 ^ (21og/3(t))A 


j=0 


j! / 


+ 


2 “ ((1-k) log/3(t)y 

2^ nt 


1-k 


j=m-l 


SO that for k>2, 


Cov {m,t,Tk) 

= am{t) - a{tf = (3{tf’"{Prn{t) “ Pitf) 


UO) _ga.2k( ik + l){2logP{t)r-^ 

g / (21og/3(t))^- / 2 \- ((l-fc)log/3(t))n \ 

/■ \1 J- J J 

In the case k = 1, we have (3'j^{t) = —(3m-i{t) with (3i{t) = 1 — t and 
/3m (0) = 1, so that 


and hence 

(50) Cov(m,t,ri) =/?(t)2(/3m(t)-/3(t)^) = -e“* ^ 

j=m+l 3 ■ 

When m = t = l, (49) and (50) yield Cov(l, l,Tfc) = —a(l)^ for all fc > 1. 
This can also be seen directly since for adjacent sites u,v in any graph, we 
have P[XY{u) = XY(v) = 1] = 0 for the annihilation process. 


4.4. The multiple annihilation proeess on T^, k >2. In this section we 
derive results for the multiple annihilation process (MAP) in analogous 
to Theorems 12 and 13 from Section 4.1, namely, (58) and (63). We consider 
only the case with k>2. 

We shall utilize the coupled construction of a realization of the MAP 
denoted {XY , t > 0) and a realization of the blocking RSA process denoted 
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> 0) in terms of an undecided site process (USP) {ZY ,t > 0) with 
state space {0,1, 2}^ , as described in Section 2.1. For any vertices v,w ^V, 
the identities (1) and (2) give us 

(51) P[zY{v)e {0,2}] =P[Y/{v)=0], 

(52) P[{ZY(v) G {0,2}} n [zY(w) G {0,2}}] = P[Y/(v) = Y/(w) = 0] 
and 


(53) P[ZY(v) G {1,2}] = P[XY (v) = 1 ], 

(54) P[{ZY(v) G {1,2}} n [zY(w) G {1,2}}] = P[xY(v) = XY(w) = 1]. 


With defined in Section 4.1, we use the following notation. For i G 
{0,1,2}, define 

■= P[zfyvo) = i], a{i,t) ■= pyj^^{v) = i\. 

The tilde indicates probabilities which refer to the USP. By (51), 

(55) /3(0, t) + /3(2, t) = Pyr^ (uo) = 0] = /3(t), 

(56) a{0,t) + a{2,t) = P[yY'"{vo) = 0] = a{t), 

where f3{t) and a{t) are as defined in Section 4.1 (and refer to the blocking 
RSA process). 

To remain “undecided” (i.e., in state 2) at time t, a site v must have > t 
and must have all of its neighbors in state 0 or 2 at time t, so that by (55), 

d(2, t) = (1 - tmo, t)+ /3(2, t)Y+^ = {i- • 


Thus by (53) and (56), we have 

(57) PlxJ^^y) = 1] = a{2,t) + a{l,t) = (1 — t)/3(t)^^^ + 1 — a{t). 


Substituting the formulae for f3(t) and a{t) derived in Section 4.1 into (57), 
we obtain the following formula for the probability that a site is occupied 
for the MAP in T^, k>2: 


(58) 


P[xY'^iv) = l] = {l-t){l + {k- i)i)(i+U/(i-fc) 


In order to calculate covariances, let v, w be vertices of with dist(u, w) = 
m. For i,j G {0,1,2}, define 


■■= P[zY^{v)=i and zf'^iw) =j]. 


For V and w to remain undecided (i.e., in state 2) at time t requires Tv > t, 
Tw > t and all neighbors in state 0 or 2 at time t. Hence by (51) and (52), 

(59) d^(2,2, t) = (1 - 


26 


M. D. PENROSE AND A. SUDBURY 


where 7 m(i),m > 0, is as defined in Section 4.1, and so is given by (37), and 
where we set 7 _i(t) = 1. 

Next, we use elementary probability followed by (51) and (52) to obtain 
amihht) = 1-P[z]^^{v) G {0,2}] -P[zf'=(u;) G {0,2}] 

(60) + P[{Z^>= (v) G {0,2}} n {Z^'^ {w) G {0,2}}] 

= 1 - 2a{t) + am{t), 

with a{t) and am{t) given in Section 4.1. 

We obtain £1^(1,2, t) indirectly, via a forward equation for 0:^(1, Ij^)- As 
h [0, we have 

P[{ZJ J^;) = Zj^^iw) = 1} \ {Zf'^ (v) = Zj- {w) = 1}] 

= 2P[{ZJ’‘ (v) = 1} n {Z^*^ (w) = 2} n {t < Tij, < t + /i}] + o(h) 

= 2h(l - tr^Pl{z[>^(v) = 1} n {Z[^ (w) = 2 }] + 0 (h), 

the last line arising because the conditional distribution of Tw, given that 
event {Z^'^ (v) = 1} fl {Z^'" (w) = 2} occurs, is uniform over the interval (t, 1]. 

Thus, 

^dm(l,l,t) = 2am(l,2,t)/(l -t). 
at 

Rearranging this identity and using (60) followed by (33) and (38), we obtain 
dm(l) 2, t) = ^ ^ ^ (—2q: (t) + a^(t)) 

(61) 

= (l-t)(/3(t)"+i-/3(t)^/3^_i(t)), 

with P(t) and (3m{t) given in Section 4.1. 

Combining (54), (59), (60) and (61), we obtain 

P[{x/’*(<.) = i}n{x?’*{u.) = i}| 

= dm(l, 1, t) + 2Q;m(l, 2, t) + 0:^(2,2, t) 

(62) . 

= 1 - 2a{t) + am(t) + 2(1 - t)/3(t)^(/3(t) - I3m-i(t)) 

+ (1 -t)^/3(t)2''7m-2(t). 

Let Covmap(w-)L^ fc) and CovR,sA(m.,t,T^) be defined as at (27), with re¬ 
spect to the MAP and to blocking RSA, respectively. In particular, CovRSA(m., t, T^) 
Oim(t) — a{t)'^. Using (62) and (57), we obtain 

CovMAp(n7.,f,rfc) 

= P[{X^Hv) = 1} n {Xf>^iw) = 1}] - (P[X^Hv) = 1])' 

= CoYKSA(m, t,Tk) + 2(1 - t)P(t)'^(a(t)/3{t) - /3m-i(t)) 

+ (1 - t)‘^P(t)‘^’'('ym-2(t) - /3(t)^). 


( 63 ) 
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In particular, Covmap(?t^, 1, = CovFi,sA(w-) 1, ?A:)j which is as it should 

be. 


5. Proof of bounds for bipartite graphs. 


Proof of Theorem 5. First we consider the blocking RSA process. 
Define the independent 1/(0,1) variables G V, by 




1 - Tv, 


if dist(u,uo) is even, 
if dist(u,uo) is odd. 


We assert that for j = 0,1,..., fc, the indicator of the event Ej, defined by 
Ej := {vj) = 0} if dist(uo,Uj) is even, and by Ej := (vj) = 1} if 
dist(uo,Uj) is odd, is an increasing function of the variables Uy,v G V. To 
see this choose any u G P, let 0 < tt < u' < 1 and fix numbers G [0,1] for 
w gV \ {u}. We assert that if 


(64) 

{Uy 

= u}n f| 

{Uw — Uw} 

implies Ej occurs, 



ui£P\{t;} 



then 





(65) 

{Uy 

= u'}n f| 

{U^ - Uyj } 

implies Ej occurs. 


u!GF\{t;} 


To see this, hrst suppose dist(u,uo) is even. Consider the effect of increasing 
the value of (and hence of r„) from u to u' while leaving the values of 
Uw,w GV \ {u}, unchanged. This can change the value of X^^ (v), if at all, 
by changing it from 1 to 0. It can also directly affect sites w neighboring v, 
if at all, by changing X^^ (w) from 0 to 1. These sites, in turn, can affect 
some of the sites tc at a distance 2 from v by changing X^^ (w) from 1 to 0. 
In general, for any w,t, the effect of changing Uy from u to u' on the value 
of xY^{w), if any, can only be to decrease it if dist(u,r(;) is even and can 
only be to increase it if dist(u,r(;) is odd. The assertion that (64) implies 
(65) follows. 

The argument is similar when dist(u,uo) is odd. In this case, increasing 
the value of Uy is equivalent to decreasing the value of Ty , and this can only 
increase the value of X^^ (w) when dist(u, w) is even and can only decrease 
the value of X^ ^ (w) when dist(u,r(;) is odd. Since dist(u,'u;) is even if and 
only if dist(uo,u’) is odd, we may conclude once again that (64) implies (65). 
This justifies our initial assertion that the indicator of event Ej is increasing 
in the variables Uy, for each j. 

The result (11) for the blocking process now follows from the Harris-FKG 
inequality. The usual version of this result (see, e.g., [10]) gives positive cor¬ 
relations of events which are increasing in a family of {0, l}-valued variables. 
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but the proof in [10] is easily adapted to the present case with [0, l]-valued 
variables Uj. 

Now consider the annihilation process. We first consider a modified ver¬ 
sion of this process, whereby for any pair of neighboring sites u,v there is 
a random time tuv uniform on [0,1] at which u annihilates v (if both sites 
were still occupied immediately before time tuv)- It is clear that only the 
event at the minimum of tuv,tvu can be effective. Therefore, this modified 
annihilation process is equivalent to the one we originally defined, except 
that for each edge the time Tg is distributed as the minimum of two inde¬ 
pendent uniform variables on (0,1), rather than as a uniform variable on 
(0,1). However, this, in turn, is equivalent to a time-change of the original 
annihilation process, as explained in Section 2.2, and therefore it suffices to 
prove the result for this modified annihilation process. 

Let the uniform random variables Uu,v be defined for each pair of neigh¬ 
boring sites u,v, by 


Uu 


tuv, if dist(u,uo) is even, 

1 — tuv, if dist(u,uo) is odd. 


By the same argument as for blocking RSA, for each j the indicator of event 
Ej is an increasing function of the variables Uuv, and the proof of (11) for 
blocking RSA carries through to the annihilation process. □ 


Proof of Theorem 6. Suppose R is a bipartite lattice and each vertex 
has degree k + 1. Let vo be some arbitrary specified vertex of V, let the 
neighbors of vq be denoted ui,..., and let the neighbors of vi (other 
than uo) be denoted rci,..., Wk- Let V* denote the graph V \ {uq}, and let 
V** denote the graph V \ {uo,ui}. 

For the blocking RSA process, define 

m-=p[xrivi)=o], 

5it) := P[xr M = xr {w 2 ) = ----.= xr* [wk) = o], 

7?(t) := P[XY^ (ui) = xY^ (U2) = • • ■ = XY^ {Vk+i) = 0]. 

If we let £ denote the event that wi,... ,'Wk are all unoccupied at time ti, 
then 5{t) is the probability of event £ for blocking RSA on R \ {ui} with the 
value of Tvq set to be (say) 2. Since Ig is increasing in 5{t) > P[£] for 
blocking RSA on R \ {ui} without any condition on Tvq. Then by positive 
correlations between sites at an even graph distance (Theorem 5), we have 

S{t)>P[£]>Pit)^ 

so that by the forward equation (Proposition 2), 

( 66 ) p>^t) = -S{t)<-P{t)Y 
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with /3(0) = 1. Hence 

d ( 1 


(67) 


dt V Pit) 


k-l / — 


<l-k 


1 


Pit)’^ 

Pit)< 


1 


Next, we assert that 


(68) r]{t)>S{t)P{t). 

To see this, note that r]it) is the probability of the intersection of events 
ACiB, where we set A:= {X^* (vk+i) = 0} and 

B := {XY^ (vi) = XY* iv 2 ) = --- = xY* ivk) = 0}, 

B' := = • • • = ivk) = 0}. 

Then P[B']=5it). 

It is to be understood that A,B,B' are all defined in terms of the same set 
of random arrival times {r^, u € H}. In effect, B' is the event that for RSA on 
V*, if we increase the value of to 2 (say) but leave r^,u G H* \ {ufc+i}, 
the same, then we have sites vi,... ,Vk vacant at time t. Then we have the 
event equality 

AnB = AnB', 


since if A occurs, then increasing the value of to 2 does not have any 
effect on the acceptances or rejections up to time t at sites other than Vk+i 
(since the arrival at site v^+i was rejected even with the original value of 
Pk+i^ at least if that value was at most t). Events A and B' are positively 
associated by Theorem 5, so (68) follows. 

Setting a{t) = P[XY (vq) = 0], by the forward equation (Proposition 2), 
along with (68) and (66), we have 

apt) = -rjit) < -Sit)P{t) = P'it)P{t). 

Since Q!(0) = /3(0) = 1, we then have 

a(t) — 1 = [ a'{s)ds< [ /3(s)/3'(s) ds = i(/3(t)^ — 1), 

Jo Jo 

so that by (67), 

a(t) < 1(1 + I3(t?) < 5 (l + (1 + ) ■ □ 
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Remark. For the annihilation process, the argument up to a point is 
the same as for the blocking RSA process. We have 


a{t) >P{t) 




m< 




2 + {k-l)t 


Unfortunately, the two inequalities appear to be in opposite directions. In 
fact, on square lattices simulations suggest that a{t) > the value on the tree 
(0.252 as opposed to 0.25). This may make this approach difficult. 


6. Proof of results relating to normal approximation. To prove Theorem 
7, we shall use the following theorem on normal approximation of sums of 
locally dependent random variables, which is a special case of Theorem 2.4 
of Chen and Shao [4]. 

Lemma 14 ([4]). Suppose V is a finite set and G V} is a collection 
of random variables, each with mean zero and finite variance. For ACV set 
Ca ■= E and suppose that Var(^\/) = 1- Aet k be a constant. 

Suppose that for each v gV there exist subsets Ay F By C Cy of V such 
that ^y is independent of fv\Av! is independent of is inde¬ 
pendent of card(C'^) < k. Then for all x G M, 

sup|P[Ck < x] - 4>(x)| < 75 k^ 
a:GM 


Proof of Theorem 7. The argument is related to that used in [1] 
and elsewhere for various geometrical central limit theorems, but here has 
no geometric content. 

Set m = card(U), and set 

r:=(51ogm with (5 :=-. 

12 log D 

For X G U, set J \ fr { v ) := N ' r { v , W ) [as dehned at (4)], and set 
St := St{V, W) = J2 (v), St := 

vGV v&V 

a := Y^Var(S't), a := \jYar{St). 

We start by observing that provided we take C{D,£) > 1, the estimate (13) 
will trivially hold if u < 1. Therefore, in the sequel we may assume with no 
loss of generality that cr > 1. 

For any two subsets U, U' of V, it is the case that J2veu and 

'fiZveU 'independent if the graph distance between U and U' 
exceeds 2r, because by definition the value of X ^^'^\ v ) is determined by 








DEPOSITION PROCESSES ON GRAPHS 


31 


{tw : w &Mr{v)} so that Suet/i® determined by the variables r^, 
w within distance r of 17, and likewise for U'. 

For V €V, set Ay '■=Af 2 riv), By '.= J\fir{v) and Cy ■.= J\f%r{v)- Also, set 

iy := - Exf^^''\v)). 

Then the conditions for Lemma 14 hold, with 


K = max(card(C^)) < (1 + Z1 + H-+ 77°'’) < D 

vSV 


6r+l 


and < d Applying Lemma 14 gives us 


^12r-+2---3, 


sup|P[(d“^(5t - ESt)) <x]- 4>(x)| < 7577°- ■ -u “m 

(69) 

< 7577^m°+^d“^. 

We need to compare St to St and cr to a. By the proof of Proposition 1, 
for any u G 17 we have 

P[X^(v) ^ xf^^'’\v)] < 777r! 

(this holds for any of blocking RSA, dimer RSA or the annihilation process). 
Hence 

(70) P[St + 7] < E[|7 -St\]< mD^/r\ < 7L[/r! 

for 7Li := e°/^77. Also, if we set Ytiv) = {X'^"'^'"\v) — X'^ (u)) we have |L)(u)| < 
1 with P\Yt{v) 0] < /r\ so that 


^Cov(Af(u),7(u;)) 

vGV w^V 


|Cov(5t,7-7)| = 

where we set K 2 := e^/^77, and similarly 
|Var(5t-5t)| = 

Hence, 


< 2m2777r! < 27L7r!, 


^ ^Cov(yi(7,yt(u;)) 

v£V w£V 


<2m^D'^ lr\<2Kllr\. 


72 - cT^I = |2Cov(5t,7 - 5t) + Var(5t - St)\< 6Kyr\, 

so that 

^ —1 <^a~^K\lr\. 

By the mean value theorem, for a G M there exists C = C{a) such that 

< max(|(l + 6 ( 7 - 27 L 27 r !)“/2 _ 17(1 _ Qa-^K^rl)^^^ - 1 |) 

< Ca-'^Kyrl 


(71) 


^1 -1 

a 










32 


M. D. PENROSE AND A. SUDBURY 


Now consider arbitrary x G M and set y := xaja. Then by the mean value 
theorem, with cf) denoting the standard normal probability density function, 
we have 


|$(y) - $(x)| <\y- x|(/>(min(|x|, \y\)) 

= (max(|x|, |y|) - min(|x|, |y|))(/>(min(|x|, |y|)) 
'a ^ a 
a 

-2 T^r 


= max 


- 1, ^ - 1) min(|x|, |y|)(^(min(|x|, | 2 /|)) 


< Ca-^Kyrl, 

where we have used (71) and the fact that max„>o{M<?i>(u)} < 1. Using (69) and (71), 
we then obtain 


|P[a-^(5i-i75i)<x]-c^(x)| 

< \P[a-\St - ESt) <y]- $(y)| + |a>(y) - ^x)\ 

< 75D‘^m^+^a-^{l + Co-’^KljE) + Ca-'^KljE. 

< (7'm^+V^, 


where in the last line we have used the assumption that cr > 1 and the (easily 
proved) fact that a <m. 

Next, for X G M we use (70) and (72) and the fact that a <m to obtain 
\P[a-\St-ESt)<x]-^{x)\ 

<\^{x + a-\ESt-ESt))-^{x)\ 

+ \P[a-\St - ESt) <x + a-\ESt - ESt)] 

-^{x + a-\ESt-ESt))\ 

< {a~^Kl/r\) + = ma~^{C'rrf + {K{/r\)a‘^/m) 

< ma~^{C'rrf + mKl/rl) = ma~^{C'rrf + K^/rl), 

since K 2 = . Hence, (72) remains valid with ESt in the left-hand side 

replaced by ESt (and also with a change of constant). Finally, by (70) and 
the fact that a <m, 

m-^a^\P[a-\St - ESt) <x]- P[a-\St - ESt) < x]| 

< m-^a^P[St + St] < ^'(KgVr!) 

(with := e^/^D) which is also bounded. Hence, having replaced ESt by 
ESt in (72), we can also replace St by St, completing the proof. □ 


Proof of Theorem 8. Let 0 < t < 1, let IT be the vertex set of a graph 
with degrees bounded by D and let U C IT be finite; set m = card(U). 


DEPOSITION PROCESSES ON GRAPHS 


33 


First consider blocking RSA. Choose an independent set I in F with at 
least m/[D + 1) elements; this is possible by a greedy algorithm. 

For V £ I, let Yy denote the indicator of the event that Ty, > t ioi w £ 
Ari(n,fF)\{n}. Then 

(73) 

VGl \ f / 

Let T be the cr-field generated by the values oi Yy,v £ I along with the 
values of 

r^„,^ceVF'\| y (A/K'y)) 

V{«e/:y„=i} 

Then, setting St := St{V,W), we have 


(74) 


YaiiSt) = YaT{E[St\T]) + .F[Var(54 |.F)] 
> F;[Var(5i|.F)]. 


If we are given the values of G /, along with the values of Tyy for w £ 
}(A/i(n))), then the remaining variability of St comes only 
from the values of Xt{v) for v £ I with Yy = l. For such n, site v contributes 
1 to St Y Ty <t and contributes 0 to St if Ty > t. Hence, 


Yax{St\Y) = Var Bin ^ T,, f = t(l - f) ^ T. 

V \t;e/ // v&I 

so that by (73) and (74), 

Var(5t) > t{l - t)Ej2yv > t{l - t)m{l - t)^/{D + 1), 

vGl 

which completes the proof of (14) for blocking RSA. 

The argument is similar for dimer RSA. Choose an independent set I' 
of edges, each incident to at least one element of V, with at least m/[2D) 
elements of I', using a greedy algorithm and the assumption that there are 
no isolated vertices in V. For each edge e £ I', let Yy be the indicator of the 
event that Tf > t for all edges / adjacent to e. The rest of the argument 
proceeds much as before. 

A similar argument applies for the annihilation process. The factor of 
(1/2) is included in the right-hand side of (15) because edges in I' might 
have only one endpoint in V. We do not give details for this case. □ 


Proof of Theorem 9. Suppose V C W, where V is finite and W is 
the finite or countable vertex set of a graph with all degrees at most D. 
Then card(A/ 5 (n, IF)) < for all v £V. Recall that IF"'' denotes the set of 
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sites in W with positive entropy. Using a greedy algorithm and the bound on 
degrees, choose a collection J of sites in U n W~^, all lying a graph distance 
at least 6 from one other, and satisfying 


(75) 


card( J) > 


card(U n W+) 


For V £ J, let Zy be the indicator of the event that for all sites w in the 
associated blocking set By (see the dehnition of positive entropy in Section 
2.7), we have Ty, < 1/2, and Tx > 1/2 for x £ (A/i(rc,lU)) \ {u;}. Then 


(76) P[Zy = l]>{lf^^^+^\ 

Let T' be the cr-held generated by the values oi Zy,v £ J, along with the 
values of 

r^,t(;GluW IJ (A7i(u)U{u}) 


Then, as at (74), setting 5i := Si{V,W) we have 
(77) Var(5i) > E[Yai{Si\P% 


If we are given the values of Zy,v £ J, along with the values of Ty, for 
u; G VF \ (U{'ugJ:Z„=i }(A/i(u) U {u})), then the remaining variability of Si 
comes only from the values of {Xy,,w £ Afi(v, W)} with v £ Jand Zy = 1. For 
such V, pick u,u' adjacent to v but not to each other. The collection of sites 
Mi{v) contributes 1 to Si if Ty < niin^^,g_A/i(i;)\{i;} Ty, [an event of probability at 
least 1/(11 +1)], and contributes at least 2 to Si if max(T„, Tu>) < minjry : y £ 
Ml (u) \ {n, u'}} [an event of probability at least 2(Z1 +1)“^]. For any random 
variable X we have Var(X) > (l/4)min(P[X > 2],P[X = 1]). Hence, 

XaT{Si\P) > i(Zl + l)-2^Z^ 

V&J 


so that by (75) and (76), 

VMS,) > 1(0 + D- T o|z„l > -dO' n 

v&J 


D^{D + iy 


which completes the proof of (16) for blocking RSA. The proof of (17) for 
the annihilation process is similar. □ 


Proof of Theorem 11. By Theorem 7, it suffices for us to prove the 
convergence of variances (19) and (20) for blocking RSA, along with the cor¬ 
responding convergence of variances for dimer RSA and for the annihilation 
process. 
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We first show the equivalence of (19) and (20). Let A" := A/'"(i) — 

Then |A”| < 1 and by the argument in the proof of Proposition 1, for any 
of blocking or dimer RSA or the annihilation process we have 


P[Ar/0]<i + 


(n + 1 — i)!' 


Since |A"| < 1, we have 


I Cov(Af, A")| < L;|A"A"| + E\/X}\E\X^\ < 2min(P[A" + 0], P[A" ^ 0]). 
Hence Var A” = 0(1) since 


n n n 

Var5^A? = j;]5;Cov(A?,A”) 

2 = 1 2 = 1 J = 1 



OO 

<32V 


2=1 
= 32e“^ 


1 


Hence, if (19) holds, then by the Cauchy-Schwarz inequality 
Cov(^^A(^,5i(L„,Z)^ =0(ni/2) 

and since St{Ln,Ln) = St{Ln,'L) +Z)r=i obtain 

VarSt(L„,L„)-Var5t(L„,Z) 

n / ^ 

(78) =Var^A” + 2Cov ^A^5t(L„,Z) 

i=l \i=l 

= 0(nV2), 

so that (19), if true, implies (20). Similarly (20), if true, implies (19). 

Define Cov(m, t) := Cov(m, f, Ti) as given by (27), with Cov(0, t) = Var Af (0). 
Then by translation invariance, 

n n 

Var(5t(Ln,Z)) = ^ ^ Cov(|z - j|, t) 

i=lj=l 
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(79) 


n—1 

= Cov(0, t) + 2n~^ ^ in — m) Cov(m, t) 

m=l 

oo 

—> Cov(0, t) + 2 ^2 Cov(m, t), 

m=l 


where the last line follows from the dominated convergence theorem. 

In the case of blocking RSA, using Theorem 13 along with the case k = l 
of (29), and collecting terms, we obtain 


OO 

^Cov(r,t) 

r=l 


(80) 


^_2t oo oo 

■EE 


i-2t) 


rH-l+2n 


2 ^ ^ (r + 1 + 2n)! 

r=l 71=0 ^ ' 

— ( E 


, even m>0 


fm\{—2t)^^ ^ /m—l\{—2t) 

odd m>0 ^ ^ 


V 2 


m: 


ml 


p-2t 

-^(—2te“^* — sinh(—2t)). 


Also, by Theorem 12, Cov(0,t) is the variance of a Bernoulli variable with 
parameter (1 + e“^*)/2, so that 


Cov(0,t) = ^(l + e-2*)(l-e-2*) 


= sinh(2t)/2, 


and combining this with (80), we see that the limit in (79) equals 
which gives us (19) as required for monomer RSA. 

In the case of dimer RSA, by using (40) and (46), and some routine algebra 
which we omit, we obtain 


OO 

Cov(0, t) +2 ^2 Cov(m, t) = 4te~^* 

m=l 


so that by (79), = 4te“^* as asserted. In fact, this can also be derived 

from the blocking RSA result by a duality argument; dimer RSA on a row 
of n vertices is equivalent to blocking RSA on a row of n — 1 vertices, where 
each site for the blocking process corresponds to a bond (edge) for the dimer 
process, and hence an occupied site in the blocking process counts as a pair of 
occupied sites in the dimer process, so that St{Ln,Ln) for the dimer process 
has the same distribution as 2St{Ln-i, Ln-i) for the blocking process, and 
hence . 

For the annihilation process, using (48) and (50) along with (79) we have 


OO 

f^AN(i) = Cov(0,t) + 2 ^ Cov(m, t) 

m=l 


oo oo 


—t 




i-ty 


m=l j=7n-\-l 
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= e-* - e-2* - 2e-* 


g g i-ty 

m=l j=m 


E 

m=l 


(-tr \ 

m\ ' 


= e-*-e-2* + 2e-*(^e-^ 

= e-2* - e-* + 2^e-2^ 

as asserted. □ 



{-ty \ 
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